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ABSTRACT 


Natural  convection  liquid  cooling  of  simulated  electronic  compo¬ 
nents  in  a  vertical  channel  was  investigated.  The  test  surface 
contained  a  single  column  of  eight  rectangular,  protruding  heated  ele¬ 
ments,  each  simulating  a  20-pin  dual-in-line  package.  Temperature 
measurements  and  flow  visualization  were  performed  for  a  number  of 
power  dissipation  levels  and  channel  widths.  Collectively,  this  infor¬ 
mation  was  used  in  interpreting  the  flow  and  transport  characteristics. 
A  correlation  to  predict  the  heat  transfer  rates  was  developed  based 
on  the  component  surface  temperatures.  Optimum  channel  widths 
were  determined  from  these  surface  temperature  measurements  for 
the  range  of  power  levels  investigated.  Temperature  distributions  in 
the  fluid  were  measured  using  a  traversing  thermocouple  probe. 
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A.  STATEMENT  OF  PROBLEM 

The  continuing  trend  toward  microminiaturization  of  electronic 
components  has  made  thermal  design  critical  to  achieving  higher 
packaging  densities.  The  ever-increasing  heat  fluxes  must  be  removed 
while  keeping  the  device  junction  temperatures  typically  below  100°C 
[Ref.  11.  The  failure  of  microelectronic  chips  increases  exponentially 
as  the  junction  temperatures  increase  [Ref.  2]. 

Microchips  act  as  discrete  sources  of  thermal  energy.  Effective 
removal  of  the  generated  energy  is  of  major  import  in  determining 
how  they  might  be  used  without  overheating  the  elements  themselves 
or  the  surrounding  package.  Thermal  control  is  one  of  the  critical 
considerations  in  the  design  of  any  electronic  cooling  system.  The 
thermal  designer’s  task  is  further  complicated  by  the  presence  of 
arrays  of  electronic  components  placed  in  proximity  to  each  other, 
requiring  detailed  considerations  of  transport  interactions  between 
neighboring  elements. 

The  design  of  packaged  electronic  components  must  consider  all 
anticipated  modes  of  heat  transfer.  A  complete  thermal  analysis  will 
account  for  all  three  mechanisms  of  thermal  energy  transport: 

•  conduction 

•  convection 


•  radiation 


The  areas  of  concern  include  transferring  heat,  via  conduction,  from 
internal  heat  sources  to  the  external  surface  of  the  package  and  the 
subsequent  removal  of  heat  from  the  external  surface  by  a  cooling 
medium.  Also,  the  coolant  temperature  must  be  maintained  at  pre¬ 
scribed  levels  for  a  given  environment. 

The  present  work  focuses  on  the  heat  removal  from  electronic 
packages  through  convective  cooling.  This  convection  may,  in  general, 
be  forced  or  natural;  the  medium  liquid  or  gas.  Because  of  the  higher 
density  and  thermal  conductivity  of  liquids  compared  to  air,  heat 
transfer  coefficients  an  order  of  magnitude  or  more  larger  than  for  air 
are  attainable.  Due  to  the  significantly  high  heat  transfer  rates  achiev¬ 
able,  cooling  of  electronic  components  by  direct  liquid  immersion  has 
been  identified  as  an  extremely  effective  thermal  control  technique 
[Ref.  2], 

Convection  heat  transfer  rate  from  a  surface  at  temperature  Ts  of 
area  A  to  surrounding  fluid  environment  T^f  is  given  by  Newton’s  law 
of  cooling; 


q'=  Ah(T,  -  rw) 

This  is  a  definition  of  h,  the  convective  heat  transfer  coefficient,  more 
than  a  phenomenological  law.  Evaluation  of  h  depends  upon  heat 
source  geometry,  orientation,  fluid  properties,  and  fluid  velocity. 

Natural  convection  results  when  fluid  density  gradients  exist  in 
the  presence  of  a  body  force  field  such  as  gravity.  Density  gradients  in 
the  fluid  commonly  occur  due  to  temperature  gradients.  Advantages  of 
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employing  natural  convection  in  electronic  cooling  include  high  relia¬ 
bility,  low  cost,  reduced  noise,  and  a  simpler  design.  Despite  these 
potential  benefits,  only  a  limited  number  of  studies  are  available  on 
liquid  immersion  natural  convection  cooling  of  electronic  components. 

Most  existing  natural  convection  studies,  both  experimental  and 
analytical,  have  considered  heating  surfaces  that  were  semi-infinite. 

Electronic  component  cooling  configurations  typically  employ  arrays  of 
small,  discrete,  flush,  or  protruding  heat  sources.  The  available  heat 
transfer  correlations  are  inadequate  in  predicting  transport  in  these 
geometries.  Convective  heat  transfer  from  an  array  of  microchips  is 
strongly  dependent  on  package  shape  and  array  density.  Indeed,  in 
real  systems,  individual  elements  or  modules  often  vary  significantly  in 
size,  shape,  and  in  the  manner  in  which  they  are  arranged  on  the 
substrate. 

The  geometric  complexity  is  further  compounded  by  the  nonuni¬ 
form  nature  of  heat  dissipation  from  the  array.  Individual  component 
dissipation  is  a  function  of  its  utilization.  Variation  in  power  dissipa¬ 
tion  from  component  to  component  on  the  same  board  and  from  board 
to  board  within  the  system  usually  exists.  .The  orientation  of  the 
boards  on  which  the  elements  are  mounted  is  yet  another  variable.  ! 

Often,  these  semi-regular  arrays  face  the  back  of  an  adjacent 
board.  The  resulting  vertical  channel  has  one  surface  which  is  rela¬ 
tively  smooth  and  the  other  one  covered  by  an  arrangement  of  large.  i 

heat-dissipating  elements.  Within  the  channel,  local  convective  heat 
transfer  is  driven  by  two  phenomena,  a  local  buoyancy  force  due  to  the 
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heated  components  and  a  forced  convection  effect  arising  from  the 
flow  history,  over  the  entire  channel  length,  including  both  heated  and 
unheated  regions.  A  simultaneous  consideration  of  these  mechanisms 
is  essential  in  obtaining  predictive  capability  for  transport  rates  in 
these  systems. 

The  complexity  of  the  problem  described  demands  an  approach 
which  builds  on  the  fundamental  physical  processes.  Efficient  cooling 
cannot  be  attained  without  understanding  the  heat  transfer  from  each 
package  by  a  determination  of  its  flow  and  thermal  fields.  The  overall 
goal  is  to  develop  a  methodology  based  on  sound  physical  principles 
which  will  allow  the  prediction  of  individual  component  temperatures, 
by  superposition,  in  a  nonuniformly  heated  array  of  components. 

B.  IMMERSION  COOLING— ANALYTICAL  AND  EXPERIMENTAL 

STUDIES 

Reference  3  provides  an  excellent  literature  review  of  previous 
work,  both  analytical  and  experimental,  in  the  field  of  liquid  immer¬ 
sion  cooling  of  electronic  equipment.  Additionally,  very  relevant  work 
has  been  performed  by  Ortega  and  Moffat  [Refs.  4-6].  They  have 
experimentally  examined  natural  convection  air  cooling  of  simulated 
electronic  components  protruding  from  a  vertical  surface,  with  and 
without  a  shrouding  wall.  Heat  transfer  rates  significantly  higher  than 
for  a  smooth  vertical  plate  were  found.  They  interpreted  this  as  an 
indication  of  turbulent  free  convection.  When  a  shroud  was  in  place, 
they  found  plate-averaged  heat  transfer  coefficients  40  to  50  percent 
higher  than  for  a  smooth,  parallel  plane  channel. 
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Jaluria  (Ref.  7]  performed  two-dimensional  computations  for 
natural  convection  air  flow  over  a  vertical  surface  with  multiple  wide, 
flush-mounted  heaters.  The  study  treated  the  flow  as  a  boundary  layer 
problem.  Rajakumar  and  Johnson  [Ref.  81  reported  computations  of 
free  convection  heat  transfer  using  rectangular  strip  heating  surfaces. 
Their  finite  element  approach  showed  significant  stagnant  areas 
between  the  blocks. 

C.  OBJECTIVE 

An  experimental  study  was  undertaken  to  to  determine  the 
hydrodynamic  and  thermal  characteristics  of  buoyancy-induced  flow 
over  several  heated  protruding  components  mounted  on  a  vertical 
surface. 

A  single  column  array  was  chosen  because  it  provided  the  sim¬ 
plest  geometry  featuring  discretized  heat  dissipation.  The  effects  of  a 
neighboring  card  were  simulated  by  a  movable  vertical  shroud  plate. 
Specific  objectives  of  this  study  were: 

•  to  visualize  the  transient  and  steady-state,  natural  convection  flow 
within  an  interrupted  channel,  for  a  range  of  component  power 
dissipation  levels  and  channel  widths; 

•  to  measure  component  temperatures  for  various  power  inputs  and 
shroud  spacings  to  develop  an  appropriate  nondimensional  cor¬ 
relation  to  predict  the  heat  transfer  rates  for  tills  geometry;  and 

•  to  measure  temperature  distributions  in  the  adjacent  fluid,  using  a 
traversing  thermocouple  probe  to  assist  in  evaluating  the  hydro- 
dynamic  and  thermal  nature  of  the  flow. 
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This  study  was  a  follow-on  to  work  described  in  Reference  3  and 
is  part  of  a  continuing  effort  to  study  natural  convection  liquid  immer¬ 
sion  cooling  of  electronic  equipment. 
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n.  EXPERIMENT 

A.  GENERAL  DESIGN  CONSIDERATIONS 

Detailed  descriptions  of  equipment  design  and  construction  are 
available  in  Reference  3.  Only  a  brief  overview  is  presented  here. 

As  seen  in  Figure  1.  the  assembly  consisted  of  a  vertical  test  sur¬ 
face  which  contained  a  single  column  of  eight  rectangular  protruding 
heated  elements.  Each  protrusion  simulated  a  20-pin  dual-in-line 
package.  The  elements  were  stainless  steel  blocks  spaced  one  inch 
apart,  center  to  center,  on  a  plexiglass  substrate.  A  12.70  mm  thick 
plexiglass  shroud,  when  used,  simulated  the  reverse  side  of  a  second 
printed  circuit  board. 

The  blocks  (Figures  2  and  3)  had  imbedded  thermocouples  on 
each  of  the  five  surfaces  in  contact  with  the  fluid  to  provide  surface 
temperature  measurements.  A  sixth  thermocouple  was  placed  in  the 
center  of  each  element  mounting  slot  for  measurements  of  heater 
temperature.  A  miniature  thermocouple  probe  mounted  on  a  three- 
dimensional  traverse  enabled  measurements  of  temperature  distribu¬ 
tions  within  the  fluid. 

A  foil  heater  mounted  on  the  component  side  facing  the  substrate 
was  powered  using  a  regulated  D.C.  power  supply.  Each  heater  was 
run  in  series  with  a  precision  resistor  and  all  eight  heaters  were  in 
parallel  with  the  power  supply.  The  current  to  each  heater  was 
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0.02  R  (4-PLACES) 
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calculated  by  subtracting  the  heater  voltage  from  the  source  voltage, 
then  dividing  by  the  resistance  of  the  precision  resistor.  Power  dis¬ 
sipation  through  each  element  was  determined  by  calculating  the 
product  of  heater  voltage  and  heater  current. 

Figure  2  shows  the  orientation  of  the  coordinate  axes  used  in  the 
discussion  to  follow.  The  front  face  of  the  component  is  that  surface 
in  the  x-z  plane  in  contact  with  the  fluid.  The  side  faces  are  the  two 
x-y  surfaces,  right  and  left  as  determined  by  viewing  the  component 
from  the  front.  The  top  and  bottom  faces  are  those  in  the  y-z  plane, 
downstream  and  upstream  of  buoyant  fluid  flow,  respectively. 
Component  height  is  in  the  x-direction,  depth  in  y,  and  breadth  in  the 
z-direction. 

The  immersion  fluid  was  purified  water  contained  in  a  one  cubic 
meter,  plate  glass  tank.  Suspended  in  the  water  were  particles  of 
Pliolite,  an  inert  pigment  with  a  specific  gravity  of  0.93.  The  particles 
were  illuminated  with  an  eight  milliwatt  helium  neon  laser  split  into  a 
plane  by  use  of  a  cylindrical  lens.  The  plane  of  visualization  was  varied 
by  realignment  of  the  laser  lens  assembly  and  the  camera. 

B.  EXPERIMENTAL  PROCEDURE 

1-  Established  Method 

Details  are  again  provided  in  Reference  3,  hereafter  referred 
to  as  the  previous  study.  Noteworthy  differences  are  noted  here. 

2.  Instrument  Settings 


Power  settings  of  0.2,  0.5,  1.0,  1.2,  and  1.5  watts  were  used 
for  each  channel  width.  Because  of  the  inability  to  de-aerate  the 


immersion  fluid,  1.5  watts  was  the  upper  limit  of  power  available. 
Settings  above  this  generated  bubbles  on  the  block  surfaces. 

3.  imtnimsBt  Bsafliaga 

When  measuring  steady-state  surface  temperatures,  succes¬ 
sive  temperature  measurements  at  10-minute  intervals  were  com¬ 
pared.  Only  when  all  thermocouples  varied  less  than  0. 10°  C  were  the 
final  measurements  taken.  The  temperature  acquisition  program  of 
the  previous  study,  modified  to  reflect  the  numbering  of  elements 
from  bottom  to  top  (flow  direction)  and  to  allow  data  storage 
(Appendix  A),  was  then  executed. 

When  using  the  traversing  thermocouple  probe,  the  data 
acquisition  system  was  programmed  to  continuously  monitor  the  volt¬ 
age  registered  by  a  single  channel.  Observation  of  the  on-screen  print¬ 
out  allowed  determination  of  temperature  variation  from  point  to  point 
within  the  flow  field,  as  well  as  at  any  point  over  time. 

A  program  (Appendix  A)  which  monitored  from  one  to  eight 
thermocouples  over  time  was  used  when  making  a  record  of  the  tran¬ 
sient  temperature  response  of  the  surface  thermocouples  to  power 
step  up  or  step  down.  Typical  data  presented  were  for  the  front  face 
temperatures  of  the  various  components. 

4.  Photographic  Technique 

Best  results  were  obtained  when  using  an  f  stop  of  2.8  and 
exposure  time  of  20  to  40  seconds.  Camera  placement  was  approxi¬ 
mately  six  inches  from  the  tank  wall  at  a  height  which  allowed  all 
eight  blocks  in  the  field  of  view.  The  camera  was  inclined  at  an  angle 
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of  about  5  degrees  from  the  test  surface  when  a  shroud  was  present. 
This  was  helpful  in  minimizing  reflected  light.  A  shroud  spacing  of 
two  times  component  height  from  the  plexiglass  substrate  was  deter¬ 
mined  to  be  the  minimum  at  which  useful  pictures  could  be  obtained. 

C.  DATA  ANALYSIS 

While  Qcond  and  Qconv  were  calculated  exactly  as  recommended 
by  the  previous  study,  further  data  reduction  differed  significantly. 

For  each  block,  a  surface -averaged  temperature  was  generated  by 
multiplying  the  temperature  and  area  of  each  fluid -exposed  face  and 
dividing  the  summation  over  the  block  by  the  total  exposed  face  area 
of  the  block.  The  average  over  the  eight  blocks  of  the  surface-averaged 
temperatures  was  then  averaged  with  the  immersion  fluid  tempera¬ 
ture  to  produce  a  film  temperature.  Fluid  properties  based  on  curve 
fits  (Appendix  B)  of  data  from  Incropera  and  Dewitt  (Ref.  9J  were  then 
calculated  for  the  film  temperature. 

A  new  characteristic  length,  L  ,  for  determining  a  nondimen- 

0 

sional  temperature  excess,  a  modified  Grashof  number,  Gr*,  and 
Nusselt  number,  Nu,  was  chosen.  The  progression  of  its  use  for  each 
block  is  as  follows: 
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The  various  symbols  used  are  defined  in  the  List  of  Symbols  section. 

Other  characteristic  lengths,  L  to  determine  block  spacing  along 
the  plexiglass  substrate  and  D  to  represent  the  horizontal  breadth  of 
the  element,  were  also  chosen.  Their  use  is  as  follows: 
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~  =  nondimensional  location  of  component  centers  along  the 
L  test  surface 

~  =  horizontal  position  along  the  element’s  front  face 

i 

Quantitative  development  was  best  expressed  as  plots  of: 

•  Nu  vs.  Gr*  relationships 

•  nondimensional  block  surface  temperature  excess  vs.  element 
position 

•  temperature  distributions  in  the  fluid  adjacent  to  the  test  surface 

•  block  surface  temperatures  vs.  time 

Pertinent  trends  in  these  variations  are  discussed  next. 
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A.  ALL  ELEMENTS  POWERED,  STEADY  STATE 
1.  Flow  VUnftHgfltlnn 

Photographs  of  the  natural  convection  flow  in  a  plane  (x-y) 
geometrically  centered  through  each  component  are  shown  in  Fig¬ 
ure  4.  The  visualizations  represent  all  elements  heated  to  power 
settings  of  0.2,  0.5,  and  1.0  watts,  left  to  right,  respectively.  The 
exposure  times  were  40,  20.  and  20  seconds,  left  to  right.  The  pho¬ 
tograph  toning  techniques  required  for  mass  printing  of  this  work  do 
not  allow  sufficient  resolution  for  detailed  study.  Should  the  reader 
desire  to  see  prints  of  much  greater  clarity,  contact  Y.  Joshi  at  the 
address  given  in  the  distribution  list  of  this  study. 

The  flow  appears  to  consist,  in  each  case,  of  strong  upward 
motion  driven  by  the  buoyant  force  near  the  test  surface.  The  fluid 
tends  to  follow  the  contours  of  the  protrusions.  The  particle  traces 
^  are  closer  passing  the  element  front  face  as  compared  to  the  inter¬ 

block  space.  Dead  zones  are  not  obvious  except  perhaps  for  very  small 
areas  where  the  blocks’  top  faces  intersects  the  substrate  and  for  a 
slightly  larger  area  above  the  top  face  of  the  last  block  in  the  array. 
There  is  no  evidence  of  recirculation  or  formation  of  dead  zones  near 
the  bottom  face  of  any  block. 

i 
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Figure  4.  Steady  Flow  In  X-Y  Plane  at  Power  Levels  of  0.2,  0.5, 
and  1.0  Watts.  Powers  Increase  From  Left  to  Right 
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In  all  three  photographs,  a  developing  boundary  layer  characteristic  is 
evident  farther  out  from  the  protrusions.  It  attains  a  thickness  of 
about  3.5  component  depths  for  0.2  watts  and  2  for  1.0  watt.  Its 
origin  moves  upstream  with  increased  power,  from  about  1  com¬ 
ponent  height  at  0.2  watts  to  3  component  heights  at  1.0  watt.  The 
velocities  appear  greater  and  the  boundary  layer  becomes  thinner  as 
power  is  increased.  Ambient  fluid  is  entrained  over  the  entire  length 
of  the  layer. 

Figure  5  presents  flow  visualization  in  the  plane  (x-z)  passing 
through  the  front  faces  of  the  elements.  All  elements  were  powered 
to  1.2  watts  and  exposure  time  was  40  seconds.  The  photograph 
shows  an  upward-moving  buoyant  layer  near  the  components.  Ambient 
fluid  is  entrained  from  each  side  into  the  flow,  resulting  in  an  increase 
in  the  horizontal  extent  of  the  upflow  region  downstream.  The  flow  in 
the  vicinity  of  the  elements  appears  three-dimensional. 

As  with  the  side  view,  the  flow  appears  to  follow  the  contours 
of  the  projections  to  about  the  same  degree  (dipping  in  about  3-5  mm 
from  the  side  faces),  the  flow  also  appears  to  wrap  around  the  cor¬ 
ners.  The  particle  traces  about  20  percent  along  the  block  (z  direc¬ 
tion)  from  the  side  edge  appear  to  be  largely  in  the  x-y  plane  and 
those  at  the  side  face  are  in  the  x-z  plane.  The  particle  traces  in 
between  are  3-D,  each  progressively  moving  along  the  90°  arc. 

Overall,  the  flow  can  best  be  described  as  following  the  con¬ 
tour  of  an  hourglass  split  in  the  x-z  plane.  The  elements  are  the  upper 
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and  lower  vessels.  This  configuration  is  repeated  over  the  length  of 
the  array.  It  becomes  more  pronounced  as  the  flow  moves  up  the 
array  and  is  eventually  shed  as  a  plume  above  the  last  element. 

2.  Quantitative 

Figures  6  and  7  are  nondimensional  representations  of  tem¬ 
perature  excess  for  0.2  and  1 .0  watts,  measured  on  various  component 
surfaces  vs.  block  position.  Also  shown  are  the  data  for  heater 
temperature  measurements.  The  two  plots  follow  the  same  patterns. 
The  numbers  for  the  0.2  watts  are,  as  expected,  larger  due  to  the 
inverse  relationship  with  T*.  which  is  h  dependent.  Also  as  antici¬ 
pated,  the  heater  excesses  are  much  greater.  The  drop  in  the  heater 
temperature  between  blocks  6  and  7  could  perhaps  result  from  a 
slightly  smaller  (~4%  less)  power  achieved  by  block  7.  Unfortunately, 
the  heater  thermocouple  for  block  8  was  inoperative.  For  both  power 
settings,  the  imbedded  surface  thermocouples  showed  increasingly 
greater  excesses  along  the  array. 

Of  primary  importance  in  thermal  design  is  the  availability  of 

the  appropriate  heat  transfer  relationships.  As  discussed  earlier,  such 

correlations  for  discrete  protruding  heat  sources  in  liquids  have  not 

been  available  in  literature  up  to  now.  To  this  end,  a  logarithmic  plot 

* 

of  Nu  vs.  Gr*  (Figure  8)  was  generated  with  eight  curves  representing 
element  position.  Each  curve  had  five  points  for  the  power  settings 
0.2,  0.5,  1.0,  1.2,  and  1.5  watts  resulting  in  a  range  of  Gr*.  Although  a 
linear  pattern  with  Gr*  was  evident,  the  spread  in  data  indicated  an 


additional  dependence  on  element  location.  A  review  of  the  data  sug¬ 
gested  a  method  by  which  it  could  be  collapsed.  A  modified  Nusselt 
number.  Nu*,  was  calculated  by  use  of  a  scaling  factor  based  on  the  net 
convected  energy  of  the  flow  at  a  downstream  location.  Since  each 
block  dissipates  approximately  an  equal  amount  of  power,  the  net 
convected  energy  is  proportional  to  the  number  of  upstream  compo¬ 
nents.  The  scaling  factor  was  of  the  form  Bna  where  Bn  is  the  block 
number  and  a  >  0.  A  value  of  a  =  1/6  was  found  to  provide  the  best 
correlation  for  the  data.  As  expected,  the  strong  three-dimensional 
effect  on  transport  results  in  a  weaker  downstream  dependence  than 
the  1/5  exponent  in  the  Fujii  and  Fuji!  correlation  for  a  uniform  flux, 
semi-infinite  vertical  surface.  Therefore, 

Nu *=  Nu  ■  Bn* 

for  each  element  at  each  power  setting.  The  modified  data  are  plotted 
in  Figure  9  along  with  the  linear  least  squares  fit.  Based  on  this,  the 
following  correlating  equation  for  the  heat  transfer  rates  is  obtained: 

Nu*=  L88(Gr  *)  '5 

for  50  <  Gr  *  <  1000  and  Pr  -  7. 


B.  ALL  ELEMENTS  POWERED,  TRANSIENT 

1.  Flow  Visualization 

Figure  10  is  a  series  of  photographs  taken  during  the  0.2  watt 
step-up  transient.  For  each  picture,  the  exposure  time  was  20  sec¬ 
onds.  Exposure  of  the  first  picture  began  with  power  application.  The 
others  were  initiated  at  20,  80,  120,  and  200  seconds,  respectively, 
following  the  start  of  heating.  Clearly,  the  first  two  show  predom¬ 
inantly  parallel  flow.  The  next  pictures  represent  increasing  devel¬ 
opment  of  the  contour  flow  and  boundary  layer  characteristic 
development.  By  200  seconds,  almost  fully  developed  flow  is 
recognized. 

Figure  1 1  is  a  similar  series,  again  of  20-second  exposures, 
taken  during  the  step-down  transient  from  1.0  watt.  The  first  visual¬ 
ization  at  far  left  began  at  power  down.  The  remaining  photographs 
were  initiated  80.  160,  240,  and  440  seconds,  respectively,  following 
shutdown.  The  steady-state  flow  characteristics  persist  to  an 
increasingly  lesser  degree  for  the  entire  period  observed.  The 
entrainment  of  ambient  fluid  lessens  at  large  times.  Additionally,  flow 
origin  moves  slightly  downstream  over  time. 

2.  Quantitative 

The  1.0-watt  step-up,  step-down  transient  temperature 
responses  for  the  component  front  surfaces  are  presented  on  the 
same  graph  (Figure  12).  Steady  state  is  achieved  about  25  percent 
faster  for  the  power  up.  The  curves  representing  blocks  1  through  8 
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excluding  3  travel  closely  together,  especially  in  the  first  quarter  of 
the  response  period  for  the  step.  The  almost  single  line  corresponds 
to  the  parallel  flow  pattern  seen  in  the  photographs.  This  indicates  a 
purely  diffusive  transport  response  for  a  large  part  of  the  transient. 
Later,  convective  effects  become  appreciable  as  entrainment  develops 
downstream.  For  the  step  down,  steady  state  temperature  is  reached 
much  more  quickly  than  breakdown  of  the  established  flow.  The  per¬ 
sistence  of  flow  is  conducive  to  the  rapid  cooling  of  the  blocks. 

The  plot  for  0.2  watts  (Figure  13)  reaches  steady  state  at 
about  the  same  time  for  both  step  up  and  step  down.  This  period  is 
longer  than  for  either  1.0-watt  case  and  indicates  that  the  stronger 
flow  more  readily  establishes  the  final  thermal  pattern. 

C.  TEMPERATURE  MEASUREMENTS  IN  THE  FLUID 

Characteristics  of  the  thermal  transport  in  the  fluid  were  mea¬ 
sured  using  a  traversing  thermocouple  probe  adjacent  to  all  protruding 
elements  and  for  points  midway  between  the  components  (Figures  14 
and  15).  The  fluid  temperature  responses  near  the  blocks  and  inter¬ 
block  spaces  follow  the  same  pattern,  with  temperature  levels  in  the 
vicinity  of  blocks  4  and  6  surprisingly  greater  than  for  the  others. 
Indeed,  both  graphs  can  be  seen  to  have  four  groupings.  Block  1  by 
itself,  the  lowest  curve,  is  followed  by  2,  3,  and  5;  then  4,  7,  and  8; 
and  finally  block  6  by  itself.  More  detailed  local  measurements  are 
needed  to  explain  these  variations. 
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The  layer  between  the  blocks  is  slightly  thicker  with  respect  to 
that  at  the  block  faces.  This  might  be  expected  from  the  flow 
visualization  (Figure  4)  showing  that  the  particle  traces  spread  out 
between  elements.  In  each  case,  the  thermal  layer  is  about  one 
component  height  or  less  and  is  therefore  deeply  imbedded  within 
the  momentum  exchange  region.  This  is  expected  due  to  the 
relatively  large  Pr  ~  7  for  water. 

Next,  the  traversing  probe  was  placed  in  the  thermal  layer,  1  mm 
from  the  surface,  and  moved  in  the  x-z  plane  across  the  test  surface 
face  (Figures  16  and  17).  The  temperature  variations,  with  position, 
reinforce  the  photographs  in  suggesting  strongly  three-dimensional 
transport.  Ambient  temperature  is  the  lowest  point  shown.  The 
peaks  at  about  15  percent  from  each  edge  (Figure  16)  correlate  with 
the  regions  in  the  visualization  (Figure  5),  where  the  flow  begins  to 
shift  from  x-y  dominant  to  x-z  and  may  indicate  a  common  region 
through  which  those  particles  are  flowing.  The  result  is  more, 
warmer  fluid  convected  over  that  region.  Likewise,  the  peaks  of  Fig¬ 
ure  17  might  represent  a  similar  flow  shift  at  locations  between  the 
components. 

Finally,  the  probe  was  placed  above  the  center  of  the  top  face  of 
block  8,  2  mm  out  from  the  substrate  to  measure  the  downstream 
temperature  decay  (Figure  18).  The  extremely  high  temperature  in 
the  immediate  vicinity  of  the  top  face  suggests  a  dead  zone,  as  did  the 
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visualizations,  where  a  small  pool  of  fluid  is  heated  primarily  by  ther¬ 
mal  conduction  from  the  component  upper  surface.  There  was  no 
sign  of  a  vortex. 

All  measurements  with  the  probe  were  remarkably  consistent 
over  time  (±1  |xv)  at  any  position.  No  evidence  by  any  measurement 
indicated  the  presence  of  turbulent  flow  for  the  power  levels 
investigated. 

D.  INDIVIDUAL  ELEMENTS  POWERED 

Individual  blocks  were  powered  to  1.0  watt,  one  at  a  time.  The 
resulting  steady-state  array  temperatures  were  monitored  and  the  flow 
visualized.  These  are  discussed  next. 

1.  Flow  Visualization 

A  series  of  photographs  show  single-element  heating  as  that 
element  changes  from  one  block  to  the  next  up  the  array.  The  expo¬ 
sure  times  for  all  pictures  presented  were  40  seconds.  Figure  19  pre¬ 
sents  left  to  right  blocks  1,  3,  and  4  heated  and  Figure  20  shows 
blocks  5,  6,  and  7.  The  flow  with  only  the  lowest  block  heated  resem¬ 
bles  the  all-components-heated  case,  indicating  flow  around  block 
contours  and  a  boundary  layer  characteristic  over  the  entire  length.  As 
progressively  higher  blocks  are  heated,  a  quiescent  area  begins  to 
develop  adjacent  to  the  lower  part  of  the  test  surface.  Typically,  fluid 
entrainment  results  in  vigorous  flow  around  the  first  block  upstream  of 
the  one  heated.  At  two  blocks  upstream,  a  comparatively  weaker  flow 
is  seen  and  at  three  the  fluid  is  nearly  quiescent. 
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2.  Quantitative 


Nondimensionalized  temperature  excesses  for  the  faces  and 
heater  of  each  element  are  displayed  (Figures  21  through  25).  The 
plots  show  very  similar  temperature  excess  patterns  around  the 
heated  element  as  the  selected  component  progresses  up  the  array. 
The  protruding  elements  have  minimal  effect  when  they  are  part  of 
the  unheated  starting  length. 

E.  COMBINATIONS  OP  ELEMENTS  HEATED 

As  in  Section  D.  progression  was  made  up  the  array,  except  that 
all  upstream  blocks  remained  heated.  Figures  26  through  30  indicate 
very  little  effect  of  powering  of  downstream  components  on  the  tem¬ 
perature  responses  of  already-powered  upstream  components.  Also 
evident  is  a  very  sharp  decay  in  temperature  in  the  unheated  blocks. 
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Figure  22.  Nimdimcnsional  Temperature  Excess  Levels  at 
Right  Faces  of  Individually  Heated  Components 
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Figure  23.  Nondimcnsional  Temperature  Excess  Levels  at 
Top  Faces  of  Individually  Heated  Components 
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Figure  24 


NO  SHROUD,  SIDE  TEMPS 


Figure  28.  Nondimcnsional  Temperature  Excess  Levels  at 
Top  Faces  of  Multiple  Heated  Components 
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NO  SHROUD,  BOTTOM  TEMPS 
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A  plexiglass  plate  (shroud)  was  used  (Figure  1)  to  simulate  the 
channel  formed  by  two  adjacent  parallel  printed  circuit  boards.  The 
test  surface  as  before  has  protruding  heated  elements.  The  shroud 
simulates  the  reverse  side  of  the  second  board  in  an  array. 

A.  ALL  ELEMENTS  POWERED.  STEADY  STATE 
1.  Flow  Visualization 

Photographs  presenting  natural  convection  flow  in  the  x-y 
plane,  geometrically  centered  through  each  component  (Figure  31), 
were  evaluated.  The  visualizations  are  for  all  elements  heated  to  0.5, 
1.0,  1.2,  and  1.5  watts.  The  exposure  times  were  20  seconds.  The 
shroud  distance  was  18  mm  from  the  substrate  or  approximately  three 
component  depths.  The  first  six  blocks  are  clearly  visible  in  each  pic¬ 
ture.  Figure  32  shows  the  top  portion  of  the  test  surface  for  0.5  watts. 
The  other  power  settings  do  not  differ  significantly  and  are  not  shown. 

The  flow  near  the  elements  appears  similar  to  that  observed 
in  the  absence  of  the  shroud.  However,  the  shrouding  wall  prevents 
ambient  fluid  from  entrainment  into  the  buoyant  up-flow  in  this  plane. 
As  power  is  increased,  the  paths  described  by  the  particle  traces  are 
more  closely  spaced,  indicative  of  increased  velocity  and  increased 
fluid  flow. 
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Figure  32.  Steady  Flow  in  X-Y  Plane.  0.5  Watts.  18  mm 
Channel  Spacing.  Top  Blocks 
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Figure  33  is  the  front  view  of  flow  in  the  x-z  plane  through 
the  front  surface  of  the  elements.  It  is  for  1.0  watt,  18  mm  shroud 
spacing,  and  exposure  time  of  40  seconds.  Unlike  the  no-shroud  con¬ 
dition,  flow  is  more  restricted  to  the  sides  of  the  protruding  blocks. 
The  characteristic  hourglass  flow  seen  for  infinite  channel  width  is 
weak  in  the  vicinity  of  the  middle  two  blocks  and  absent  elsewhere. 
The  entrained  fluid  turns  very  sharply  into  the  flow  direction.  Because 
the  flow  is  restricted  to  the  channel  width  and  is  predominant  on  the 
block  sides,  its  velocity  is  much  greater  here.  Consequently,  there  is 
more  forceful  entrainment  of  fluid  in  the  x-z  plane. 

2.  Quantitative 

Due  to  the  large  amount  of  data  collected,  temperature 
excesses  for  power  settings  of  1.2  and  0.2  watts  only  are  presented  in 
graphical  form.  These  settings  are  representative  of  the  upper  and 
lower  ends  of  the  power  range  available.  Channel  widths  presented 
are  6.2,  8,  14,  18,  and  23  mm  measured  from  the  substrate. 

The  first  width  was  just  barely  greater  than  component 
height.  This  spacing  is  representative  of  pure  conduction  dominated 
transport.  The  nondimensionalized  temperature  excess  vs.  element 
position  graphs  for  0.2  watts  (Figures  34-38)  and  1.2  watts  (Figures 
39-43)  show  significant  difference  between  the  smallest  spacing  and 
all  others.  With  that  spacing  excluded.  Figures  44  through  53  present 
the  relevant  data  in  a  manner  better  suited  for  analysis.  The  plots  for 
0.2  watts  show  the  next  closest  spacing,  8  mm  from  the  substrate,  to 


55 


58 


59 


62 


64 


66 


^  -o 


>C  +  < 


•^CE  < 


0.2  0.4  0.6  0.(1 


X4-  3C 


73 


<  + 


XD  <  +; 


S  -Kl 


Q 

-cs3:2w^ 


x  c  « 


x  HSl 


a  cix 


SHROUD 


6 


yield  higher  temperatures  over  the  range  of  element  positions.  For 
higher  spacings,  there  was  no  consistency  in  the  relative  position  of 
the  curves  from  one  element  position  to  the  next.  This  indicates 
relative  independence  of  transport  on  channel  width  for  spacings 
greater  than  8  mm.  However,  data  for  1.2  watts  includes  even  the 
8  mm  spacing  within  the  mix  of  curves  (its  curve  does  not  consistently 
yield  higher  temperatures).  Therefore,  for  1.2  watts  the  indepen¬ 
dence  of  transport  extends  below  the  8  mm  spacing.  General  patterns 
observed  for  both  powers  are: 

•  a  steady  rise  in  front  and  side  face  temperatures  with  element 
position  as  expected  due  to  exposure  to  increasing  convected 
energy  at  higher  positions. 

•  a  relative  levelling  of  top  and  bottom  temperatures  at  higher  posi¬ 
tions.  This  is  perhaps  due  to  increasing  fluid  velocities  along  the 
channel. 

•  the  drop  in  heater  temperatures  at  block  7,  which  has  four  per¬ 
cent  less  power  on  average. 

The  differences  between  the  two  powers  result  from  the 
stronger  driving  force  at  1.2  watts.  Increased  fluid  velocities  result  in 
improved  heat  transfer  characteristics,  making  8  mm  an  acceptable 
spacing.  This  point  is  emphasized  by  referring  to  the  graphs  of  com¬ 
ponent  face  temperature  excess  vs.  shroud  spacing  for  0.2  watts 
(Figures  54-58)  and  1.2  watts  (Figures  59-63).  Significantly,  for  the 
1.2  watts  there  is  no  further  temperature  reduction  for  a  spacing 
greater  than  8  mm.  Indeed,  there  is  a  slight  rise  in  temperature  of 
the  latter  blocks  (7  and  8)  when  the  shroud  is  removed.  These  higher 
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temperatures  (at  the  last  blocks)  are  not  observed  when  a  shroud  is  in 
position  because  the  global  forced  convection  within  the  channel  is 
driving  the  fluid  more  strongly  past  the  last  block.  As  seen  from  Fig¬ 
ures  54  through  58,  face  temperatures  at  a  power  level  of  0.2  watts  for 
8  mm  are  somewhat  higher  than  for  greater  spacings.  By  14  mm,  the 
difference  with  larger  spacings  begins  to  become  indistinguishable. 
Clearly  for  this  geometry  and  for  the  power  settings  examined,  a 
channel  width  of  2.5  times  component  depth  is  most  satisfactory 
while  a  width  of  1.3  times  element  depth  may  be  acceptable.  Further 
confirmation  of  the  observation  that  spacings  beyond  8  mm  do  not  sig¬ 
nificantly  affect  heat  transfer  is  seen  in  Figures  64  through  68.  These 
show  the  Nu*  vs.  Gr*  data  along  with  the  correlation,  Nu*  =  1.88 
(Gr*)15,  obtained  in  Chapter  III.  The  correlation  for  an  unshrouded 
surface  holds  very  well  through  14  mm.  At  8  mm  there  is  a  slightly 
greater  deviation.  As  anticipated,  it  is  inapplicable  for  the  smallest 
spacing  of  6.2  mm. 

B.  INDIVIDUAL  ELEMENTS  POWERED 

As  in  the  absence  of  a  shroud,  individual  blocks  were  powered  to 
1.0  watt,  one  at  a  time,  but  with  a  shroud  spaced  14  mm  from  the  sub¬ 
strate.  The  resultant  data  is  displayed  in  Figures  69  through  73.  The 
plots  follow  similar  trends  as  for  the  no-shroud  condition  (Figures  21- 
25)  but  show  somewhat  smaller  temperatures.  Again,  the  protruding 
elements  have  minimal  effect  on  the  unheated  starting  length. 


!  64.  Nu*  vs.  Gr*,  With  Correlation.  23  mm  Channel  Spi 
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C.  COMBINATIONS  OP  ELEMENTS  HEATED 

In  a  manner  identical  to  that  without  a  shrouding  wall,  progres¬ 
sion  of  elements  heated  was  made  up  the  array  while  all  upstream 
blocks  remained  powered.  Additionally,  block  8  was  powered 
throughout  and  a  shroud  14  mm  from  the  substrate  was  in  place.  The 
plots  (Figures  74-78)  reveal  that  for  this  channel  width,  as  with  no 
shroud,  the  powering  of  downstream  components  has  minimal  effect 
on  those  upstream  components  already  powered.  The  plots  further 
show  a  weak  effect  by  a  heated  element  (block  8)  on  those  compo¬ 
nents  between  it  and  the  ones  powered  upstream,  which  are  heated 
only  by  the  convective  flow  generated  upstream  of  them. 

D.  ALL  ELEMENTS  POWERED,  TRANSIENT 
1.  Flow  Visualization 

Figure  79  shows  transient  flow  development  in  the  x-y  plane 
following  1.0  watt  power  step-up  with  a  19  mm  spaced  shroud  in 
place.  All  exposures  are  for  20  seconds.  The  photographs  from  left  to 
right  are  initiated  at  0,  20,  40,  100,  and  140  seconds  following  start  of 
heating.  In  the  first  picture,  parallel  flow  is  developing  over  blocks  3 
through  8  with  relatively  little  movement  at  blocks  1  and  2.  By  picture 
number  3,  the  flow  follows  more  closely  the  contours  of  the  protrud¬ 
ing  blocks  and  originates  further  upstream.  By  photograph  5,  the  flow 
appears  nearly  fully  developed. 

The  power  step-down  from  0.2  watts  was  again  for  a  19  mm 
channel  width  and  20-second  exposures  (Figure  80).  The  pictures 
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Figure  78. 
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from  left  to  right  were  initiated  at  0,  60,  200,  320,  and  460  seconds, 
respectively,  following  cessation  of  power.  As  with  infinite  channel 
width,  the  flow  characteristic  persists  over  the  entire  period  observed. 
Only  in  the  last  picture  does  the  flow  begin  to  follow  less  closely  the 
contour  of  blocks  7  and  8  and  originate  farther  downstream. 

2.  guantiiativfi 

The  transient  temperature  response  curves  for  powers  of  1.0 
and  0.2  watts  and  a  shroud  spacing  of  19  mm  (Figures  81  and  82}  fol¬ 
low  exactly  the  same  pattern  as  those  for  no  shroud  (Figures  12  and 
13).  For  the  spacing  examined,  the  shroud  does  not  affect  the  trans¬ 
port  evolution. 
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Figure  81.  Transient  Temperature  Response 
Following  Initiation  of  Heating.  1.0  Watt 
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Figure  82.  Transient  Temperature  Response 
Eolfowlng  Cessation  of  Heating.  0.2  Watts 
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V. 


Flow  visualization  in  concert  with  quantitative  analysis  provided  a 
number  of  significant  determinations: 

•  The  convective  flow  was  three-dimensional. 

•  The  flow  was  laminar. 

•  The  thermal  layer  was  imbedded  within  the  momentum  exchange 
region. 

•  Heated  elements  have  minimal  effect  on  those  elements, 
upstream,  which  are  already  heated. 

•  There  is  no  gain  in  heat  transfer  for  channel  widths  greater  than 
two  times  component  depth. 

•  A  correlation  has  been  found  to  predict  heat  transfer  rates  for  the 
subject  geometry. 
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VI. 


It  is  suggested  that  the  following  areas  of  study  be  experimentally 
explored: 

•  Using  the  same  test  surface: 

-  explore  transient  response  for  other  block  faces 

-  with  the  traversing  thermocouple  probe,  determine  variation  of 
temperature  with  position  (in  the  thermal  layer)  for  all  blocks 

•  Construct  new  test  surface  and  rerun  the  same  battery  of  experi¬ 
ments  to  determine  significant  changes  when: 

-  heated  elements  are  flush  mounted 

-  heated  elements  are  spaced  closer  in  the  x-direction. 
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SOFTWARE 

TEMPERATURE  ACQUISITION  PROGRAM 
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PRINT  " 

BLOCK  83" 

630 

PRINT  "  “ 

640 

OUTPUT  709 (“ CONFHEftS 

OCV. 210-215. USE  0“ 

650 

FOR  1-30  TO  35 

660 

ENTER  709i Volts! I ) 

670 

Tenp! I )-. 0006797+!  25825.  !328*Volts< I >  >-< 607789.2467*!  Volts!  I )• Volts! I ) > ) 

2  1 952034. 3364* <  Volts!  1 1*3) >♦! 8370810996. 1874*! Volts! I >“4>  > 

680  PRINT  MT.  C.  r'tl+l,"  Volts  O.C.  ” i Volts! I > ."Tenp.  DE6.  C  "(Tenp!  I > 

690  NEXT  I 

700  FOR  >1  TO  14 

710  PRINT  •*  “ 

720  NEXT  J 

730  PRINT  "  BLOCK  8Z" 

740  PRINT  "  " 

750  OUTPUT  709 \ "CONFHEftS  OCV , 2 1 6-2 1 9 , USE  0" 

760  FOR  1-36  TO  39 

770  ENTER  709(Volts!I> 

780  Tenp! 1 >-.0006797+!  25825.  13Z8*Volts! I ) )-! 607789 . 2467-! Volts! I )♦ Volts! I > )  )-! 
2 1952034. 3364*!  Volts!  I ) *  3) )+!  8370810996. 1 874* ( Volts! I >’4)> 

790  PRINT  “T.  C.  8"(I+1."  Volts  O.C.  ” (Volts! I ) ."Tenp.  DEG.  C  " i Tenp! I ) 

800  NEXT  I 

810  OUTPUT  709 i "CONFHEftS  DC V .300-301 .USE  0“ 

820  FOR  1-40  TO  41 

830  ENTER  709 (Volts!  I) 

840  Tenp!  I )-. 0006797+!  25825. l328»Voltsi  I > ) -! 607789 . 2467*! Volts! I >*Volts! I  )  )  )-! 
21952034.3364*!  Volts!  I ) *  3) )  +  < 8370810996. 1874*! Volts! I  )*4) ) 

850  PRINT  "T.  C.  8"|I+1."  Volts  D.C.  “  ;  Volts!  I  >  ."Tenp.  0E6.  C  "  i  Tenp!  I  > 

860  NEXT  I 

870  PRINT  "  " 

880  PRINT  '■  BLOCK  81" 

890  PRINT  "  “ 

900  OUTPUT  709 ( "CONFHEftS  OCV . 302- 307 . USE  O" 

910  FOR  I -42  TO  47 

920  ENTER  709 (Volts!  I ) 


9  30  Te«p< I )-. 0006797+125825.  1328»Volts<  I ) )-< 607789 . 2467+1 Vol ts( D+Voltsl I  ) )  >-( 

Z 1 952034. 3384»( Voltst I >“3>  >+1 8370810996. !874»<  Voltsl I > * 4 )  ) 

940  PRINT  "T.  C.  «“il*1.“  Volts  D.C.  “ s Voltsl I > .'Tenp.  DEG.  C  "tTenplI) 

950  NEXT  I 
960  ! 

970  i 

980  OUTPUT  0Path3i Te*p< • > 

990  I 

1000  t 

1010  PRINT  "  " 

1020  PRINT  “  “ 

1030  PRINT  "  BATH  TEMPERATURES  (TOP  TO  BOTTOM)" 

1040  PRINT  "  " 

1050  OUTPUT  709 i "CONE HERS  OCV . 317-31 9 .USE  0" 

1060  FOR  I -57  TO  59 
1070  ENTER  709iVolta<I> 

1080  Te«p< I )». 0006797+1 2S825. l328»Volts( I > )-( 607789. 2467+1 Volts< I >• Voltsl I )) )-( 
21952034. 3364+ <  Volt s< I )*3) )♦( 8370810996. 1 874+ ( Volt s< I )*4>  > 

1090  PRINT  "VOLTS  D.C.  ” j Voltsl I ) .“TEMP.  OEG.  C  "»Te»p<I> 

1100  NEXT  I 
1110  END 


1  14 


DATA  REDUCTION  PROGRAM 


10  ||ll|!{!!l|l!!||i)!!|l|ll!!l!!{!l!l!l||l!!!!!!!!MMi;t!l!|imii| 

20  '  i  " 

30  m  DATA  REDUCING  PROGRAM  (BASIC)  M 

40  i  i  !  i 

50  l  »  l  i  l  i  l  l  n  i  It  M  l  i  l  l  i  l  i  H  M  M  M  M  •  M  M  M  M  M  •  M  H  H  I  M  H  M  M  H  H  l  M 

60  REAL  Tenp(  48) , Tt( 48  > , Qcond(  8) .Qconv< 8> ,Qf lux(  8 ) ,H( 8 > , Nu< 8  > ,Nustar<  8  ) 

70  REAL  Gr<  8> ,Ttt(8).At(8).Dtf( 8) ,Otr<  8).Dtl(8),Dtt<8> ,Otb(  8) ,Oth<  8 ) ,B1( 8) 

80  REAL  J*tar<  8) . Jp( 8 ) .Gratar<  8 > 

90  REAL  Power  .Space, Oxpg.Oxr  .Kpg.Kr  .Block 

100  REAL  A.G.Oia.Try ,KhZo,Beta,Kirwisc,Pr 

110  INTEGER  I.J.K.L.M.N.O 

120  1 

130  Oxpq-. 006731 
140  Dxr-.003175 

150  Kpg-.1421 

160  Kr“ .0389 
170  Chpht-. 00787 
180  Chpdpth" . 006 1 0 

190  Chpudth-,02388 

200  Pr ie«<  2»(  Chpht+Chpudth) )♦<  2«< 2«< Chpht+Chpdpth) > >♦<  2»(  2*<  Chpudth+Chpdpth) ) ) 

210  A1  •Chpud'th*Chpht 

2Z0  AZ"Chpud-th«Chpdpth 

230  A3«Chpht*Chpdpth 

240  Atotal“A1 ♦<  2 . »AZ  >♦<  2 . »A3) 

250  Atotala“<  2 . »A2 )  +  (  2 . *A3) 

260  G-9.81 

270  |  M  t  I  I  !  I  !  !  !  !  1  M  M  M  M  !  M  M  H  M  M  M  M  M  I  M  M  M  I  1  '!  M  •  M  I  M  I  I  I  I  It  M  M  •  M  M 

280  M  CALCULATE  CHARACTERISTIC  LENGTH  " 

230  i  i  i  i  ;  i  i  ;  i  it  i  i  i  i  |  ii  i  i  i  i  i  m  1  '  *  1  i  i  i  i  i-t  m  i  i  i  i  i  l  i  i  i  i  I  I  i  m  i  M  l  i  ii  i  i  i  i  I  |  i  l  M  M  i 

300  Lbar»Atotal/Prin 
310  ' 

320  R1»Oxpg/<Kpo«AI  ) 

330  RZ-Oxr/(Kr«A1  ) 

340  i 

350  PRINTER  IS  1 

360  PRINT  "INPUT  RUN  #,  INTEGER  ONLY" 

370  INPUT  Try 

380  PRINT  “  " 

390  PRINT  "  " 

400  PRINT  "INPUT  POWER  VALUE  IN  WATTS" 

410  INPUT  Power 

420  PRINT  "  ” 

430  PRINT  "  " 

440  PRINT  "INPUT  BLOCK  #’S  POWERED" 

450  INPUT  8Inp$ 

460  PRINT  "  “ 

470  PRINT  "  " 

480  PRINT  "INPUT  SPACING  IN  MM" 

490  INPUT  Space* 

500  PRINT  "  " 
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SI®  PRINT  "  “ 

52®  PRINT  "IWHJT  AM8IENT  TEMPERATURE" 

S3®  INPUT  Tinf 

S4«  PRINT  “  11 

SS«  PRINT  “  11 

S60  «  PRINT  "INPUT  NEW  FILENAME*1* 

S7®  !  INPUT  Neufile* 

58®  (  PRINT  “  11 

59®  *  PRINT  ”  11 

E®®  PRINT  “INPUT  OLO  FILENAME*11 

61®  INPUT  Oldfile* 

620  I 

63®  I 

640  !  CREATE  BOAT  Neufile*. IZ8 

6S0  I  ASSIGN  ORoad  TO  Neufile* 

660  ! 

670  ASSIGN  BStreet  TO  Oldfile* 

680  ENTER  •Street  * Te*p< » ) 

690  FOR  I-®  TO  47 

700  J-J+1 

71®  (PRINT  J.Te«p<I> 

72®  NEXT  I 
730  ! 

74®  ! 

750  ( l ( l ( l ( ( I | f | ( I t ( ( M ( I i l I  I l | I  I | I ( I  I l ( l t  l  I  (  |  M (  I  ( ( ( ( M  (  l  (  l  l  l  I  I  (  l  M  I 

760  » »  CALCULATE  AVERA6E  SURFACE  TEMPERATURES  ! • 

770  I  |  l  j  |  |  I  I  |  I  (  I  |  |  M  1  I  1  M  1  M  !  H  M  I  !  1  M  I  M  |  (  |  l  I  )  1  M  1  (  (  1  M  I  (  (  1  *  1  H  M  !  !  ( 

78®  !  PRINT  "AV6.  SURFACE  TEMPERATURE11 

79®  J-0 

800  FOR  1-1  TO  8 

81®  IF  1-6  THEN 

82®  At< I )•( < Tenp<  J+ 1 )»A3  >♦<  Te*p(J>2 )*A3>+( Tenp(  J+3  >+AZ )+(  Tenp<  J+4 )*AZ  )  \! 

Atotala 

830  ELSE 

84®  At< I )-< ( Tenp< J1+AI )  +  <  Teap<  J+l >»A3>+( Te*p( J+Z > »A3>+< Tempi  J  +  3>«AZ  >♦<  Te 

mp<  J+4 )»AZ ) )/Atotal 
850  ENO  IF 

86®  i  PRINT  At<I) 

870  J-J+6 

880  NEXT  I 

89®  i 

90®  II  II  II  M  I  I  I  I  I  I  II  I  I  I  I  I  I  I  I  I  I  I  I  II  I  II  |  II  I  II  I  |  II  H  I  I  I  II  II  I  It  II 

91®  M  CALCULATE  TFILM 

9Z0  I l | M l i I l i l i l i M  |  l i t i i l II i i i i I i l l I | i  |  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 
93®  J-l 

940  Tbar»<  At( J)+At( J+1 )+At(  J+Z )+At(  J+3)  +  At( J+4)+At<  J+5>+At( J+6)+At<  J+7)  >/8.® 
950  Tf il«c-< < Tbar+Tinf )/2.®) 

96®  Tf ilm-Tf ilnc+273. IS 
97®  i  PRINT  “TFILM  DEGREES  C-.Tfiimc 

98®  !  PRINT  “TFILM11  ,Tf  ila 

990  i  i  1 1 1  |  i  i  1 1  n  n  1 1 1  i  i  i i i 1 1 ii i 1 1 1 1 1 1 1 1 1 1  |  i  |  1 1 1 1  i  |  i  m  i|  i  i  |  i  m  |  i  m  i  i  i  i 

1®®®  ! !  CALCULATE  UATER  PROPERTIES  <  ! 


i  i  i  i  i  i  i  i 
i  i 

i  I  i  i  i  i  i  i 
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1010  !  M  I  l  !  !  !  <  M  M  I  •  M  f  M  !  !  M  M  •  M  M  *  •  t  !  !  II  •  •  l(  !  M  M  I  •  II  M  I  M  !  I  M  i  I!  •  ! 

1020  8eta-<-7.9448548E-8*<  Tf il«*2 > >♦< 5. 73S64?9E-S*Tf il*>- . 0097810SE3 
1030  Spvol-(  4.69639E-9»<  Tfiln*2>>*<-2. S374S£-6*Tf iln>+. 00) 341903 
1040  Oynvxsc-< 3. 234851 IE-7«< Tf  tla‘2 > >♦< - .0002 U74487»Tf xl*>+.  036166792 
1050  Kxnvlsc-Spvol*Oynvisc 
1060  Kh2o-(  1.4l818l818E-3«Tf xl*>+. 1866 

1070  Pr-( 4.65706E- 3*< Tf il«*2>>*< -2.922094»Tf il»>+463.  3319 

1080  ! ! ! ! !  !  i  l  If • ! ! ! ! ! ! ! ! 1! H ! ! ! ! ! ! 11 ! ! ! 1 ! ! • M ! !  !  !  ! !  !  I M ! • ! ■ ! I l |  |  I  |  |  l  l  I 


1090  M  CALCULATE  TENP-TINF  M 

1  100  !!l!!!!!!!l!!lll!!!l!!l!l)lll!!l!!!!!l!l!!!l!!!!!!!!!!!l!!l!!!!!| 

1110  !  PRINT  “ TEMP- T INF" 

1120  FOR  I«0  TO  47 

1130  IF  1-5  THEN 

1140  Tt(  I >-888.888 

1150  ELSE 

1160  IF  1-30  THEN 

1170  Tt(  I >-888.888 

1180  ELSE 

1190  Tt< I )-Te«p( I )-Tinf 

1200  END  IF 

1210  ENO  IF 


1220  •  PRINT  Tt( I ) 

1230  NEXT  I 

1 240  !!!!!!i|!!l!!lll!!ll!!!MM{!!M!!!!l!!i!!!!!H!!M!!l!!!!MII!!M 

1250  1!  CALCULATE  CONOUCTION  LOSSES  THR0U6H  THE  TEST  SURFACE  !! 

1260  !  M  (I  l  M  M  M  •!  I  I!  M  M  M  M  M  M  M  H  M  M  l  !  M  I  M  H  M  M  M  M  M  M  M  1  M  !  M 

1270  !  PRINT  “QCOND" 

1280  L*S 

1290  FOR  1-1  TO  8 
1300  IF  1-1  THEN 

1310  Qcond< I >- . 015»Pouer  *- 

1320  ELSE 

1330  Qcondt I )«Tt<L)»< 1.0/IR1+R2)) 

1340  END  IF 

1350  i  PRINT  Ocond(I) 

1360  L-L+6 

1 370  NEXT  I 

1  380  I  l  l  M  M  M  M  M  l  M  M  !  !!  I  H  M  '  '  M  I  I  •  M  I  *  i  M  •  i  i  l  I  •  I  l  M  M  •  !  M  I  I  M  I  I  M  M  I 

1390  »'  CALCULATE  CONVECTED  HEAT  FLUX  <i 

1 400  l l | I  I i i i | i t l l i M ! ! I  I  I  I  *  M I  I ! I! ! I  I  I ! I  1 M ! 1 H ! I i  !  I  I  I  I  I  I  I  I  I  I  I M 1  I • I ! I  I 
1410  i  PRINT  "QCONV" 

1420  FOR  1-1  TO  8 

1430  Oconv< I )-Pouei — Qcond< I  ) 

1440  Qf lux(  I )>Qconv( I > / A 1 
1450  NEXT  I 

1 460  I  I  I  I  !  l  l  l  M  |  ! M ! i • i • ! i I  I M I  * i I  I i * i I i ! • i i ' 1 ! I  I  *  * ! i • ! I i I  1 1 1 ! I M  I  I  •  1  I  M 

1470  H  CALCULATE  TEMPERATURE  SCALING  FACTOR  * * 

1 480  !  l  |  !  '  !  I  !  I  I  1  •  M  M  I  !  I  i  *  I  '  I  !  i  !  I  M  M  I  1  I  M  l  M  I  1  M  I  !  !  i  *  M  '  I  I  I  i  *  M  •  l  I  •  M  l  ! 

1490  I  PRINT  "T*" 

1500  FOR  1-1  TO  8 


1510  T#tar<I>-Oflux<I>»Lbar/KhZo 
1520  !  PRINT  TSTAR(I) 

1530  NEXT  I 

1 540  II M  t i | I  I i l | | l| ' M  t  t  M I ! • l l l l | I  II i  i  l  l  It  I  l  M  |  M  ll  M  t l !  (  i  i  i  M  i  |  (I  •  i  l  i 

1550  !  !  CALCULATE  MODIFIED  GRASHOF  NUMBER  I ( 

1 560  l  It  (I  •  !  II  !  !  H  !  I  !  I  H  M  I  I  !  !  •!  H  H  !  1  1  •!  !  !  !  !  !  !  ■  !  !  H  !  M  M  !  !  M  I  t  i  i  H  I  !  !  »  * 

1570  FOR  1-1  TO  8 

(580  6rstar< I )-G»Beta«Qflux< I >•< Lbar*4 >/(  Kh2o»<  K inviac*2  > > 

1590  NEXT  I 

1600  l|  I  I ! H  I  I  I  I  I  1 ! I  I H I  I  I  I  I  I t I  I  '  I  I  I  I  I  M  I  H  I  I  I  I  H  1  I!  ! I  I  I  I  1  I  I  I  I i i  I  i  M  I  I  I  I 

1610  H  CALCULATE  HEAT  TRANSFER  COEFFICIENT  II 

1 620  |  |  |  I  |  II  I  I  I  I  I  I  I  I  I! I!  I  1  I  I  !  M  I  1  I  I  I  I  !  I  '  I  '  I  I  I  I  H  I  M  !  I  1!  I  I  I  I  I  I  I  I i I  I  I  I  I  I  1  I 

1630  FOR  1-1  TO  8 
1640  H< I )-Qf lux(  I )/<  At(  I )-Tinf ) 

1650  NEXT  I 

1 660  !!!!ll!!HIHI!!!!U!l!!l!)!l!ll!((!(!l!(!l!!(!l!!!!!ll!IMii!M!!!! 

1670  II  CALCULATE  NUSSELT  NUMBER  (1 

1 680  !IH(!MIIII!!M!!!(!l!!!(!!l!ll!IMIH!!lll!ll!(ini!mui!ll!H|! 

1690  FOP  I-l  TO  8 

1700  Nu(  I )-H< I )»Lbar/KhZo 

1710  NEXT  I 

1720  I  I  I ! I  I ! ! ! I  I ! ! 11 1 ! I ! I  I ! II  I i I ! I  I  I  I M  I  *  I  I  !  I  M  I  H  M  !  I  I  I  !  I  I  I  I  l *  1!  !  I  ) M  M  I 

1730  I  I  CALCULATE  NUSSELT  SCALING  FACTOR  I  I 

1 740  l I  II  I  I M !! I  I • I  II ! I  I  I  I ! I  I M I  I  I M  M •  I  I  I  I  •  l  I  I  I  I  I  !  I  1  I  I  I  I l l  I  I •  •  l  I  I  l  I  I  l  l  I  I 
1750  FOR  1-1  TO  8 
1760  Js-9-I 

1770  Nuatarl  I >-Nu< I )»Js*< 1/6) 

1780  NEXT  I 

1 790  M  I  I  I  I  II II 1 l I  I  I  I  I  II  I M I  I  I M I  I  i  M  1  '  M  *  I  I  M  I  I  I  H  I  II  1  M  '  I  I  •'  i  i  II  i  i  ll  M 

1800  I!  6ENERATE  OUTPUT  DATA  FILES  'I 

1810  I  I  I  l  III  I  |  |  Ml  I  |  ll  I  ill  1 1  M  I  I  I  1  l  M  '•  I  I  I  I  I  1  I  I  »  M  I  l  I  I  I  M  |  |  i  l  l  l  |  l  M  l  I  l  l  l 

1820  I 

1830  PRINTER  IS  701 

1840  I 

1850  PRINT  USING  . RUN  NUH8ER"" , 1 SX ,00" ; Try 

1860  PRINT  USING  . POUER  IN  WATTS"" , 1 1 X , Z. D” (Power 

1870  PRINT  USING  . SPACING  IN  MM"" , 1 ZX , I0A" s Space* 

1880  PRINT  USING  . BLOCK  #’s  POWERED"". BX, ?0A"iBlnpS 

1890  PRINT  USING  . AMBIENT  TEMP.  OEG.  C"" ,5X,DD. DO" i Tinf 


1900  PRINT  USING  . FILM  TEMPERATURE"" .9X.D0D. DO" » Tfiln 

1910  PRINT  USING  . FILM  TEMPERATURE.  C”" , bX ,00. DD" »  Tt line 

1920  PRINT  USING  . THERMAL  COND.  OF  WATER"" . 3X . Z . DOD" i KhZo 

1930  PRINT  USING  . EXPANSION  COEFFICIENT , B"" . ZX . 306" t Beta 

1940  PRINT  USING  . KINEMATIC  VISCOSITY"" ,6X. 30E" iKinvisc 

1950  PRINT  USING  . PRANOTL  NUMBER"" . 1 IX ,00.00"iPr 

I960  PRINT  "  " 

1970  PRINT  "  " 

1980  PRINT  "  " 


1990  PRINT  "BLOCK  #" ."OCONO" ."QCONV" ,"GR  #" ,"T*" ,"H" ,"NU  *","NU» 
2000  FOR  1-1  TO  8 


2010  Block-9-I 

2020  PRINT  USING  "3X .0, 3X , S3DE , 3X, SZ .000, 3X ,S30E  . ZX ,S3DE , 3X , S3DE . 2X , S30E , 2X  .S3 
DE" i Block , Ocond( I ) ,  Qconvl I ) .Srstarl I ) .Tatar! I > . H< I ) ,Nu< I ) .Nustarl I ) 

2030  NEXT  I 

2040  PRINT  ”  M 

2060  PRINT  “  “ 

2060  PRINT  "  " 

2070  PRINT  “  “ 

2080  !!!l!!ll!!!!!!!llM|!!!!|||!!!!!!!l!!!!l!lMl|l!!!!!!<!!!i|!!!lMII|! 


2090  I •  CALCULATE  TEMPERATURE  DIFFERENCES  ! I 

2100  !  !  |  |  |  !  |  !  *  I  I  1  !  !  !  ! ! ! ! ! I M ! 1 ! ! ! I ! I  I ! M ! ! < ! ! II ! M ! • ! ! M  f ! *  I ! M  M ! <! ! I  I ! M 

2110  K-0 

2120  FOR  N- I  TO  8 
2130  Dtf(Nl-Tt(KI 

2140  OtrINl-TtIK-H  1 

2150  Dtl( N)-Tt( K+2 ) 

2160  Ott<N)-Tt<K+3> 

2170  Dtb<Nl-Tt<K+4> 

2180  Oth<N1-Tt<K+5i 

2190  K-K+6 

2200  NEXT  N 

2210  PRINT  "BLOCK  #“  ,“OTF"  ,“OTR"  ,“DTL''  .“DTT"  ."DTB"  ,“DTH" 

2220  FOR  1-1  TO  8 
2230  Block -9-1 

2240  PRINT  USING  “3X,0, 3X.S30E , 3X, S30E , 3X , S3DE , 3X ,S3DE . ZX.S3DE , ZX , S3DE" {Block , 


Dtf < I ) ,  Dtr( 1 1  ,Dt  1<  1 1  ,Dt t<  1 1  ,Otb<  1 1 ,Dth< I > 

2250  NEXT  1 

2260  !!!|!!l!MMII!!l|li!!!l|!|l|l|ll'!||l||llMllill!lH!IMiliil|imH 

2270  ! t  CALCULATE  NON  DIMENSIONAL  TEMPERATURE  DIFFERENCES  ' ! 

2280  I  l  M  M  I  M  M  l  M  II  I  M  M  M  M  M  •  M  !  M  l  M  !  M  !  t  •  M  M  |l  »  •  II  !  M  l  II  l  i  M  l  l  ll  ll 

2290  FOR  I»1  TO  8 

2300  !  IF  1-1  THEN 

2310  '  Othnl I  1-88888 

2320  '  ELSE 

2330  '  IF  1-6  THEN 

2340  I  Difnl I  1-88888 

2350  Dtfnl  1 1  -D-tf  <  I  l/Tstarl  1 1 

2360  Otrn(  I )-Otr< I l/Tatarl  1 1 

2370  Dtln<  I  )»Otl(  I  >/T»tar<  I  > 

2380  Ottn( I l«Ott( I l/Tskarl  I ) 

2390  Otbn( I !-Otb( I I/Tatar! 1 1 

2400  Othn< I )«Oth< 1 1/Tstar(  1 1 

2410  IF  1-1  THEN 

2420  Othnl  1 1-88888 

2430  ELSE 

2440  IF  1*6  THEN 

2450  Dtf n< I  1-88888 

2460  ENO  IF 

2470  END  IF 

2480  NEXT  I 

2490  PRINT  "  “ 

ZS00  PRINT  "  " 


1  19 


ZSI0  PRINT  “X/L"  .  "OTFN"  .  "OTRN"  .  “OTLN"  ."DTTN"  .  “OTBN" .  “OTHN" 

2520  Xs-,795 

2530  Sap* .0254/. 3048 
2540  FOR  I-t  TO  8 
2550  X1-X9-((I-1 >»6ap> 

2560  PRINT  USIN6  ''.DOO,3X  3DE , 3X , S3DE , 3X , S3DE , 3X , S3DE ,2X , S3DE , 2X , S3DE" j  XI ,Dtf 
n< I ) .Otrn< I ) ,Otln< I > ,Ottn< I .  Hbn<  I )  .Dthn<I) 

2570  NEXT  I 
2580  PRINT  "  " 

2590  PRINT  “  “ 

2S00  PRINT  “OTR  -  RIGHT  T.  C.  -  TINF.  etc." 

2610  PRINT  “T»  IS  TEW>£RATUR£  SCALIN6  FACTOR" 

2620  PRINT  "SR  *  IS  NOOIFIEO  GRASHOF  NUMBER- 
2630  PRINT  "H  IS  HEAT  TRANSFER  COEFFICIENT- 
2640  PRINT  "NU  8  IS  NUSSELT  NUMBER" 

2650  >  OUTPUT  SRoadi Qcond< • > , Oconv(  • ) ,Gr< • ) , Ttt< • > ,B1(  • ) ,Dtf  (  » ) ,Dtr( • ) , 

Otl( « ) ,Ott< • ) .Dtb< » ) .Dth(  * ) 

2660  * 

2670  I i  NOTE  THIS  PROGRAM  IS  USED  IN  CONJUNCTION  UITH  THE 

2680  H  THE  DATA  AQUISITION  PROGRAM  AND  ACCOUNTS  FOR  THE 

2690  ! !  INOPERATIVE  THERMOCOUPLES  ON  BLOCKS  3&6 

2700  ! 

2710  END 
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THERMOCOUPLE  VARIATION  WITH  TIME  PROGRAM 


10  |  t  H  l  t  H  l  I  I  II  !  I  M  H  '  M  |  M  l  |  I  I  U  l  M  ■  M  l  l  M  I  M  l  |  M  H  M  i  <  M  M  i  M  M  l 

20  ! !  TEMPERATURE  TRANSIENT  RESPONSE  PR06RAH  ■ » 

30  | ! li i | i | l ii ! | |  !  i  l  i  |  |  !  |  I  I  i  l  |  I  !  I  I  i  M  l  |  I  l  |  |  II  !  I  l  l  |  |  !  I  l  I  l  it  i  i  |  |  i  i  l  i  | 

*0  REAL  T2(  575) , T3<  575 ) , T4( 575 ) , TS( S75) , T6<  575) ,T7(  57S ) , T8( S7S1 
S0  REAL  V2<  575 )  , V3( 575  > ,  V4<  S75) ,V5( 57S> ,  V6( 575 ) ,  V7<  S75 ) ,  V8<  575) 

G0  REAL  6< S7S ) 

70  N-S7S 
80  PRINTER  IS  701 

90  CREATE  BOAT  "TRANS _ :  ,700,0, 1"  .502 

100  ASSIGN  0Trai  11  TO  "TRANS _ 700,0,1“ 

1  |0  m  i  i  l  l  l  |  |  |  i  i  l  l  ill  i  l  M  t  l  l  l  m  i  i  i  l  t  i  i  i  t  i  n  i  i  i  i  l  i  i  i  l  l  i  |  l  i  |  i  i  |  ii  i  i  ii  l  i  i  m  i  m  i  i  i 

120  BEEP 

130  FOR  J-l  TO  5 

140  OUTPUT  709 ;  "CONFHEAS  DCV _ USE  0“ 

1S0  ENTER  709 1 V 1 ( J ) 

1G0  OUTPUT  709 1  "CONFHEAS  DCV _ USE  0" 

170  ENTER  709 i V2( J) 

180  OUTPUT  709 1  "CONFHEAS  OCV _ USE  0" 

190  ENTER  709»V3( J) 

200  OUTPUT  709 i" CONFHEAS  DCV _ USE  0" 

210  ENTER  709 1 V4<  J ) 

220  OUTPUT  709 1  “CONFHEAS  OCV _ USE  0“ 

230  ENTER  709jV5<J) 

240  OUTPUT  709 j" CONFHEAS  DCV, _ USE  0" 

250  ENTER  709»VG(J) 

260  OUTPUT  709 i "CONFHEAS  OCV _ ,USE  0" 

270  ENTER  709|V7<J)  *~ 

280  OUTPUT  709i  "CONFHEAS  OCV _ USE  0" 

290  ENTER  709|V8<J) 

300  NEXT  J 

310  BEEP 

320  FOR  J-G  TO  N 

330  OUTPUT  709  i  "CONFHEAS  OCV _ USE  0" 

340  ENTER  709 i V 1 ( J ) 

350  OUTPUT  709  i  "CONFHEAS  DCV _ USE  0" 

3G0  ENTER  709iV2<J) 

370  OUTPUT  709 t "CONFHEAS  DCV _ .USE  0" 

380  ENTER  709|V3(J) 

390  OUTPUT  709 1 "CONFHEAS  DCV, _ USE  0" 

400  ENTER  709,V4(J) 

410  OUTPUT  709  i  “CONFHEAS  DCV _ USE  0" 

420  ENTER  709tVS<J> 

430  OUTPUT  709 1  "CONFHEAS  DCV _ USE  0" 

440  ENTER  709|VS(J) 

450  OUTPUT  709 1" CONFHEAS  OCV _ USE  0" 

4G0  ENTER  709 i V7<  J  > 

470  OUTPUT  709 1" CONFHEAS  DCV _ USE  0" 

480  ENTER  709sV8<J) 

490  IF  J-N  THEN 

500  BEEP  2450.97»Tone, 1 .5 

510  END  IF 

520  NEXT  J 
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530  A-. 0006797 

540  8-25825.  I3Z8 

550  C— 607789.2467 
560  D— 21952034.  3364 
570  E -83708 10996. 1874 
580  FOR  J-l  TO  N 

590  IF  J<5  THEN 

600  6<J)-0 

610  ELSE 

620  6<  J)-< i-5)»T in 

630  !  Tin  IS  TIME  BASED  ON  TIME  PROGRAM  RUNS  FROM  2ND  BEEP  TO  TONE  DIVOED 

640  !  BY  N-5.  E.6.  FOR  7  THERMOCOUPLES  Tin  IS  .745,  FOR  8  IT  IS  .852 

650  ENO  IF 
660  NEXT  J 

670  PRINT  "  " 

680  PRINT  n  “ 

690  FOR  J-l  TO  N 

700  T1<  J)-A+(B»V1<  J)  )*<C*<  VI  <  J)*Zm<0*<  V1(  J)*3>  )+(£*<  VI  <  J)*4>> 

710  T2<  J  )-A+<  B»W2<  J )  )+<  C»<  V2<  J )  *2  ) )+(  D«(  W2(  J  >“  3> )  +  (  E*<  VZ(  J )  *4 )  ) 

720  T3<  J)-A+<B»V3(  J)  )+(C*(  V3<  J  )*2  >  )♦<  D»<  W3(  J  >  “  3 »  )♦(£*<  V3<  J)*4>  > 

730  T4< J)-A+(B»V4( J)  >+(C»(W4< J) *2 ) >♦(  D«< V4< J)*3) )  +  <£•( V4( J) *4) > 

740  TS< J)-A+(B«V5< J))+< C*(WS( J)*Z ))♦<□•(  V5<  Jr 3) )♦<  E»( V5<  J) “4 ) ) 

750  T6<  J)-A-KB«V6<  J)  )  +  <  C«(  W6<  J  )“2  )>+(D*(  V6(  J  >*3>  >♦<£•<  V6(J)*4>  > 

760  T7< J)«A*(B*V7< J))+(C»( W7( J ) *2 ) >+< D»< V7(  J ) “ 3 > >+< EM V7< J)*4) ) 

770  T8< J)-AMB«V8<  J >)MCMV8<  J )*2 ) )+< D«< W8(  J)’ 3) >+<£•< V8(  J>*4> ) 

780  PRINT  USING  ” 1 X .DOO . 1 X .4DE . 1 X ,40E .  1 X ,40E . 1 X . 40E . 1 X . 4DE , 1 X ,4DE . 1 X , 4DE . 1 X . 4 
OE" ;G<  J) , T 1 <  J > , T2 <  J ) ,T3<  J) ,T4<  J) ,T5( J) ,T6<  J ) ,T7(  J) ,T8( J) 

790  NEXT  J 

800  OUTPUT  BTrai 1 1 16(  •  ) ,  T1  ( •  ) .  T2< • ) . T3( • ) . T4< • ) . T5< • > . T6<  » > , T7< • ) , T8(  * > 

810  ENO 
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APPENDIX  B 


EQUATIONS  FOR  DETERMINING  FLUID  PROPERTIES 

p  =  -7.945E-10  T2FaM  +  5.736E-5  T fRM  -  9.781E-3 

v  =  4.697 E- 9  T7faM  -  2.537E-6  T FUM  +  1.342E-3 
H  =  3.235E-7  T2FaM  -  Z147E-3  Trnjt  +  0.03617 

v  ~  V/l 

k  =  1.42E-3TfflJtf  +0.187 

Pr  =  4.66E-3  TFaM  -  2.92  TFaM  +  463 

Curve  fits  were  determined  for  range  of  temperatures  from  280  to 
300  K. 
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APPENDIX  C 


SAMPLE  CALCULATIONS 


The  calculations  were  based  on  temperatures  derived  from  a  run 
with  all  blocks  powered  to  1.0  watt,  no  shroud  in  place.  When 
appropriate,  values  are  calculated  for  block  number  6. 


A.  CHARACTERISTIC  DIMENSIONS 


Perimeter  = 


2(0.00787  +  0.02388)  +  (2(  2(0.00787  +  0.00610))) 
+  (2(2(0.02388  +  0.00610))) 


=  (12393  m 


| 

A,  =  ( Q.02388X  a  00787)  =  1879  x  10"4m2 


A2  =  (0.02388X0.00610)  =  L457  x  10'V2 


A3  =  (0.00787K 0.00610)  =  4.80  x  10  5 m2 

Atotal=  0-879  x  10"4)  +  (2)  (l.457  X  10'4)  +  (2)(4.80x  10'5)  =  5.75  x  10'm2 

l 

A  tot  ala  =  (5.75  x  10‘4  -  1 879  x  10'4)  =  3.87  x  10'4m2 


r  -  5.75x  10'4 
"  0.2393 


=  240  x  10"3m 


I 


124 


B.  SURFACE  AVERAGED  TEMPERATURES 


(19.82)(  1.879  x  10"4)  +  (19.48)U  80  x  10"5) 

+  (19.49)(4.80  x  10"5)  +  ( 20.62) (L457  x  10"4) 
+  (20.26)(l457  x  10'4) 


AVG  (*) 


(5.75  x  10-4) 


-  =  20.08  °C 


(19.4 1)(  1.879  x  10~4)  +  (18.99)(  4. 80  x  10"5) 

+  (19.05)( 4. 80  x  10‘5)  +  (19.88)(L457  x  10'4) 
+  (20.08)(L457x  10-4) 


AVG  (7) 


(5.75  x  10'4) 


=  19.64  °C 


( 19. 16)(  1.879  x  10“4)  +(18.67)(4.80x  10"5) 

+  (18.6 1)(4. 80  x  10"5)  +  ( 20.00) (L457  x  10"4) 
+  (19.57)(L457x  10‘4) 


AVG  (6) 


(5.75  X  10~4) 


=  19.40  °C 


(1&96)(  1.879  x  10"4)  +  (18.33)(4.80x  10-5) 

+  (l 8.47)( 4. 80  x  10'5)+  (19.68)(L457  x  10"4) 
+  (19.37)(l457x  10"4) 


AVG  (5  ) 


(5.75  X  10"4) 


=  19.16  ^ 


( ia65)(  1.879  x  10"4)  +(18.30)(4.80x  10‘5) 

+  (18.39)(4.80x  10'5)+  (19.75)(L457  x  10'4) 
+  (19.06)(l457x  10“4) 


AVG  (4) 


(5.75  x  10'4) 


=  18.99  'X: 
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(17.76)(4.80x  10's)  +  (18.10)(4.80x  10'5) 

L+  (19.18)(l457  x  10'4)  +  (18.87) ( 1457  x  10~4)J 


'  AVG  (3) 


(  3.87  x  10"4) 


=  18.77  °C 


( 18.07)(  1.879  x  10"4)  +  (17.40)(  4. 80  x  10'5) 

+  (17.79)(4.80x  10'5)  +  (18.85)(l457  x  10'4) 
L+(18.56)(L457x  10'4) 


1  AVG  (2) 


(5.75  x  10‘4) 


=  1 8.32  °C 


(17.73)(  1.879  x  10"4)  +(17.25)(4.80x  10“5) 

+  (17.25)( 4. 80  x  10~5)  +  (18.60)(L457  x  10'4) 
+  (17.78)(  L  457  x  10’4) 


AVG  (1) 


(5.75  X  10'4) 


=  17.89  °C 


C.  TFILM 

7  =  20.08+  19.64+  19.40+  19.16  +  1 8.99  +  1 8.77  +  18.32  +  17.89  _  19  (BoC 


FILM 


19.03+  13.23 
2 


Id  13°C 


D.  TEMPERATURE  EXCESSES 

Using  data  from  block  6,  temperature  excess  is  evaluated  as: 


ATf  =  19.16-  13.23=  5.93^0 


A7,  =  18.67-  13.23=  5.44^ 
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ATt  =  20.00  -  13.23=  6.77  °C 


A Ta  =  19.57  -  13. 23  =  6.34°  C 


A7„  =24.20-  13.23=  1Q97°C 


E.  CONVECTED  HEAT  FLUX 


<7 


// 


0.984 

L  879  x  10‘4 


=  5240^— 
m 


F.  TEMPERATURE  SCALING  FACTOR 

(5240)(240  x  1Q~3) 

(Q  597) 

• 

G.  MODIFIED  GRASHOF  NUMBER 

^  (9. 8  0 ( 1 63  x  10~6)(  5240)  (  240  x  10 
(0  597)(ll2  x  10"*) 


H.  HEAT  TRANSFER  COEFFICIENT 


I.  NUSSELT  NUMBER 


Nu 


(849)(2.40x  1Q~3) 
(0.597) 


=  3.41 
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